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Abstract 

The quark mass dependences of light element binding energies and nuclear scattering lengths are 
derived using chiral perturbation theory in combination with non-perturbative methods. In par- 
ticular, we present new, improved values for the quark mass dependence of meson resonances that 
enter the nuclear force. A detailed analysis of the theoretical uncertainties arising in this deter- 
mination is presented. As an application we derive from a comparison of observed and calculated 
primordial deuterium and helium abundances a stringent limit on the variation of the light quark 
mass, 5m q /m q = 0.02 ± 0.04. Inclusion of the neutron lifetime modification under the assumption 
of a variation of the Higgs vacuum expectation value that translates into changing quark, electron, 
and weak gauge boson masses, leads to a stronger limit, \5m q /m q \ < 0.009. 
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I. INTRODUCTION 



The Standard Model is widely believed to be a low-energy manifestation of a more general 
theory that unifies the four fundamental forces of nature. Several candidate unified theories 
suggest that spatial and temporal variation of fundamental constants is a possibility, or 
even a necessity, in an expanding Universe (see, e.g. the reviews . Studies of Big Bang 
nucleosynthesis (BBN) provide a unique probe of the values of fundamental constants in the 



pre- recombination Universe. A further motivation to consider the response o 
changes in m q , the light quark mass 1 , is related to anthropic considerations 



ight nuclei to 
, |4| that have 



e.g. been used in the context of carbon production in hot stars js- 7] in order to understand 
how much fine-tuning is necessary amongst the fundamental parameters of the Standard 
Model in order to allow life to emerge on earth. Only now - based on methods as used 
here - is one really able to study the explicit quark mass dependence of the nuclear forces 
and nuclear properties and therefore their impact on, e.g., nuclear abundances, because such 
issues can only be investigated systematically and completely based on chiral effective field 
theories or lattice simulations (or combinations thereof). 

In addition, recent studies of quasar absorption spectra suggest a cosmological gradient 
in the value of the fine-structure constant, a, across the universe 8|. The existence of this 
spatial variation could be confirmed from complementary astrophysical studies such as Big 
Bang nucleosynthesis |q|. If the values of fundamental constants were different in different 
regions of space at the time of nucleosynthesis, this could be seen in the spatial distribution 
of primordial deuterium abundances. Note that while BBN is relatively insensitive to co- 
variation 10(|, the limits placed on quark mass variation in this work can be related to 
the variation of a under a range of unification models llj. Indeed many of these grand 
unification theories predict that relative variations in the strong force would be an order of 
magnitude or two larger than those of the electroweak forces (for a simple explanation of 
this see, e.g., 12]). This is also connected to anthropic questions, for if a spatial variation of 
fundamental constants were to exist, we should not ask how finely tuned the fundamental 



1 Throughout most of this paper, we work in the isospin limit m u — rrid and only consider the average 
light quark mass, m q — (m u + m^)/2; in addition the quark masses are studied in units of a fixed Aqce>, 
for only the variation of dimensionless quantities is meaningful. However, in Sec. I VIII we address the 
constraints from neutron /3-decay which necessarily requires the inclusion of isospin violation by strong 
and electromagnetic effects. 
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parameters are, but instead conclude that life emerged in the region of the universe where 
the parameters allowed for it. 

Relating the observed primordial abundances to the values of fundamental constants at 
the time of Big Bang nucleosynthesis requires theoretical models for how the nuclear reaction 
rates depend on observable quantities such as binding energies and scattering lengths, as 
well as a model for how those quantities in turn depend on the fundamental constants. The 



former has been dealt with previously in several works, see e.g. |10l.ll3Nl5| and the references 
therein; in this paper we provide a response matrix based on the method described in [3] 
for some of these quantities which are of importance to the current work. The second part 
of the problem, relating bulk nuclear quantities to values of fundamental constants, is the 
focus of this paper. 

Most of the previous studies in this context were performed on the basis of model- 
dependent estimates for quark-mass dependences of nuclear properties [16l42l|. However, 
there are two theoretical tools available that allow, in principle, for a model-independent 
access to quark mass dependences. On the one hand there is lattice QCD, on the other 
hand one has chiral perturbation theory (ChPT). The former is a direct evaluation of QCD 
in Euclidean space-time and thus the quark mass dependence is explicit. In the latter case, 
the operator structure of quark mass terms is fixed by the QCD symmetries; in fact, ChPT 
is a faithful representation of the spontaneous and explicit chiral symmetry breaking of 
QCD 22] . The strength parameters (usually called low-energy constants, LECs) of those 
operators have to be fixed from other sources. Generally, this is done by comparison with 
experiment. However, for operators that explicitely involve the quark mass, as it is the case 
here, such a determination is difficult since in nature the quark masses take definite values. 
To determine the LECs of such operators, one can either fit to lattice data directly (see, 
e.g., [23|, where the formalism is outlined for the p meson) or from low-energy scattering 
data when using some unitarization scheme in addition to ChPT (see, e.g., Refs. 24M26|). It 
should be stressed, however, that in the latter case some model-dependence is involved, since 
the quark mass terms of higher order than what was put in from ChPT are not complete and 

n 

depend on the scheme used [27[. In some cases, this model-dependence can be controlled, 
to some extent, by a comparison with lattice data. 



In Ref. 



an effective field theory treatment of the impact of quark mass variation on 



BBN was presented for the first time. In this work the quark mass dependence of the NN 
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scattering lengths was used as primary input. To constrain these, the results of Ref. 29 1 



were used, since at present the lattice is not sufficiently accurate to provide precise values 



of these fundamental parameters. Central to the analysis of Ref. 29] was a naturalness 
assumption for the quark mass corrections to the leading quark mass independent contact 
interactions. On the other hand, the same LEC was allowed to vary in a different range in 



Ref. 



• 3oi, 



which led to quite different quark mass variations of the two-nucleon properties. 



We remark, however, that the considerations in Ref. 



29l | were consistent with the earlier 



resonance saturation study of the leading and next-to-leading order contact interactions 



performed in Ref. 



3l| 



In this work we study systematically the impact of quark mass variations on two-nucleon 
observables based on a study of the quark mass dependences of mesons, since those are 
expected to give the most prominent contributions. In particular, if a strong quark mass 
dependence of, say, the potential part of the nucleon-nucleon (NN) interaction that comes 



from a exchange were present — e.g. in Ref. 16] a striking strange quark mass dependence 



of the a is conjectured — it might, in the effective field theory where this field is integrated 
out, lead to an unnaturally enhanced LEC accompanying some contact interaction. An 
example for such a pattern are some low-energy constants of dimension two that appear in 
ttN scattering: when extracted from data in standard ChPT they appear to be unnaturally 
large, however, this can be understood phenomenologically by observing that they are mainly 
dominated by the exchange of A isobars [3^]. Consequently, once the A contribution is 
subtracted the residual LECs get reduced significantly. Analogously one should expect that, 
once all meson exchange induced large quark mass effects are treated explicitely, the bulk 
of the quark mass dependence is included 2 . Still, such a procedure involves some modeling 
that induces some systematic uncertainty which is very difficult to specify. 

Our main focus here are the sigma and the rho meson. Both appear as resonances in the 
two-pion system. The cleanest way to connect their properties to the NN sector is via a 



dispersion integral of the Omnes type as used, e.g., in Refs. [35|, [36J. Here, however, we use 



a method which is technically easier to handle and more transparent, although admittedly 



In addition, we need to assume that there are no large quark masses dependences coming from sources other 



than meson exchanges. In this sense the findings of Ref. [33| are im por tant for here it is demonstrated that 
potentially large quark mass dependences induced by ttNN cuts [34j cancel, once final state interaction 
effects are considered explicitly. 
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of lower theoretical rigour: in Ref. [3l| it was shown that the four nucleon operators of the 
NN potential can be understood quantitatively in terms of the exchange of heavy mesons in 
the sense of a resonance saturation. In that paper explicit expressions are presented for this 
kind of matching. Thus we here use the following strategy: we determine the quark mass 
dependence of the light resonances using the methods of Refs. 



24 



261 ] which allows us to 



predict the quark mass dependencies of the four-nucleon contact terms using the expressions 



of Ref. 31J. To complete this study we then quantify the impact of the determined quark 
mass dependences of the mesons together with that of the nucleon, which is already studied 
on the lattice, on the NN observables via an explict calculation of scattering lengths. It is 
important to note that, to our knowledge, no explicit calculations for a dynamic generation 
of the omega- meson exist. Thus, we assume throughout that its quark mass dependence is 
the same as of the p. Clearly, this should be refined in future studies. 

The paper is structured as follows. In Section [TT] we derive the quark-mass dependence of 
nucleon and meson masses, which we use to calculate the impact of quark-mass variation on 
the two-nucleon potential in Section IIII1 The theoretical uncertainties of our calculation are 
discussed in Section HVl From the two-nucleon observables we are able to calculate the quark- 
mass dependence of helium nuclei (Section [V}. Finally, in Sections [VTl and IVHl we calculate 
the dependence of primordial abundances on nuclear observables and combine this with the 
results of the previous sections to derive a limit on the variation of the light quark-mass at 
the time of big bang nucleosynthesis. 

II. QUARK MASS DEPENDENCE OF HADRON MASSES 

Here, we study the quark mass dependence of the pertinent hadron masses. The results 
for each hadron H are most appropriately presented in terms of the dimensionless parameters 
Kh defined via 

SMh f Mr 

°-f^ = , (2-1) 

orrif rrif 

where rrif denotes the mass of the quark of flavor / and M# denotes the mass of hadron H. 
In what follows we will choose / = q for the light quarks (in the isospin limit) and / = s for 
the strange quark. Note, although rrif by itself is not renormalization group invariant, the 
quantity of relevance here, namely 8 m / /rrif is, for the same reason as quark mass ratios are 
well defined. 
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A. Quark mass dependence of the nucleon mass 



Due to significant advances in lattice QCD the pion mass dependence of especially the 
nucleon is now known to some precision. E.g. in Ref. {37I the dependence of the nucleon 
mass on the pion mass squared as calculated by the BMW collaboration is given. It is 
straightforward to extract from this the quantity K q N (for the definition see Eq. ( 12. ip ). One 
finds K q N = 0.04. Note that more recent evaluations from other lattice collaborations give 
similar results, as nicely reviewed in {ssj. It is also pointed out in that reference that the 
nucleon mass can be well represented by a linear function of the pion mass in most lattice 
simulations, which is at odds with the chiral constraints on this observable. 



Alternatively one may determine K q N from the pion-nucleon sigma term, ct^n- Ref. 



39| 



finds 



On the other hand in Ref. 



a wN = 44.91H MeV . 



(2.2) 



40 



41] a value 



a wN = 59 ± 7 MeV 



(2.3) 



is found. In what follows we use the first value as it is consistent with the classical determi- 
nation o:: 



Ref. [42[ based on dispersion relations. A completely reliable upgrade of the value 



from 42|] can only be obtained based on the recently proposed Roy-Steiner equations for 



pion-nucleon scattering that allow for precise determination of the pion-nuc 
amplitude in the physical region as well as inside the Mandelstam triangle 



eon scattering 



43 



44). There 



is also a large spread of values determined from lattice QCD which encompasses the range 



of values given above; see Ref. [45] and the recent review by Kronfeld 



46] 



Using the Feynman-Hellman theorem er/ = irifdM^ / 'dmj, one finds straightforwardly 
K N = gj/Mm and with that, based on the numbers given above, 



K 1 _ n n 4 c+0.002 



(2.4) 



consistent with the number quoted above within 2a. The values given are significantly lower 
than those presented in Refs. Q] due to the unusually large value of the tcN sigma term 



in those works. 
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B. Quark mass dependence of meson masses 



Due to their nature as (pseudo)-Nambu-Goldstone bosons (NGBs) of the approximate 
chiral symmetry of the strong interactions, the quark mass dependence of the members of 
the pseudoscalar octet is given by standard ChPT, which is model-independent. At tree 
level, for the pion, one finds 

Kl = 0.5 ; K' = 0. (2.5) 

In our calculation, we have included higher order terms in the light (u,d ) q uark mass de- 
pendence of the pion, using the SU(2) ChPT expansion up to two- loops |47j ]. The strange 
quark mass dependence of the pion and the quark mass dependence of the kaon and eta 
(which will be needed later 



br the calculation of Kjj and K[) have been calculated using 
SU(3) ChPT to one loop 48J. We remark, however, that these masses are possibly affected 
by large higher order corrections - this is an open issue in three-flavor ChPT. For our study, 
we do not want to enter this discussion but rather use the NLO corrections based on the 
standard scenario for the strange quark condensate. To get at least some feeling for the cor- 
responding theoretical errors, we have estimated the size of the higher order corrections by 
treating the ChPT expansions in two different ways: in the first of them the expansions are 
written in terms of the physical masses and decay constants of the NGB; in the second, they 
are written as a function of the tree level constants M 0n , M x and F . In our calculations, 
these tree level constants are obtained by fitting the ChPT expansions for the masses and 
decay constants of the NGB to their physical values. Since the difference between the two 
treatments corresponds to higher orders in the expansion, it serves as an estimate for the 
systematic error due to the truncation of the ChPT series. 

The values that we provide in the first line of Table [J for K( are an average of the 
results obtained using SU(2) and 577(3), the two methods of truncation mentioned above 



and different determinations for the ChPT LECs 



49 



521 ] . The error is taken so that it covers 



all the results with their statistical uncertainties, which arise from the errors of the LECs. 



Following a similar procedure we provide, also in Table HJ the KL values for the pion decay 



constant F n . 

For the other light resonances the situation is less systematic, since all of them are 
unstable and their correct description requires a pole in the complex energy plane that 
cannot be obtained within the ChPT expansion, which, up to some logarithmic terms, 
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K q 


0.494;™ 




0.00 ±0.05 


K 


0.058 ± 0.002 


K P 


0.02 ± 0.04 




0.081 ± 0.007 


K 


0.01 ±0.05 


K k 


0.048 ±0.012 


K k 


0.00 ±0.06 



TABLE I: Estimates for the K^ R coefficients and their uncertainties. 



corresponds to an expansion in powers of the energy or the meson masses. By construction, 
the amplitudes of ChPT are only perturbatively unitary and valid only near threshold. In 
Ref. [23| a formula to be used in chiral extrapolations for vector meson masses was presented, 



t 



re quark 



53, 



54 



56j, 



however, the quality of lattice data was not sufficient to pin down the slope o 
mass dependence, which is required here. Of course, there are better data now 
so one could refresh the analysis of Ref. 



23J . Here, we follow another path, which can also be 
used to explore the quark mass dependence of the rho and the sigma. Employing dispersion 
relations for the inverse mr scattering amplitude and using ChPT to fix the subtraction 
constants — where the subtraction points can be chosen in a regime where ChPT is valid - 
solves both problems and, by generating poles, allows for the study of resonances without 
a priori assumptions about their existence or nature. This is done in a way consistent with 
the ChPT expansion without introducing spurious parameters where an uncontrolled quark 
mass dependence may appear 3 . This method is called the inverse amplitude method (IAM) 
[57I ] and has been used to study both the o and the p in an SU(2) one loop treatment in 
Ref. [24] and to two loops in {25 1. The SU(3) version of this study can be found in Ref. 26]. 

Let us note that, within the IAM, all the dependence on the quark masses appears 
through the NGB masses, which are explicitly present in the amplitudes, both kinematically 
and in interaction vertices. Thus, we can calculate the K* R parameters for the resonances 
generated within the IAM by varying the masses of the NGBs, whose dependence on the 
quark masses was discussed in the previous section, and measuring the corresponding change 
on the position of the poles. 

We have performed this calculation using SU(2) and SU(3) and different sets of LECs 



3 Clearly, the expansion is controlled only to the order of the chiral expansion used as input. Terms of 
higher order produced by the formalism are not necessarily correct [27(, although, at least they will 
respect unitarity and the correct analytic structure of the amplitude. 
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obtained from I AM fits 24-26 



In each case, we have changed the masses of the 
NGB using the two different methods for the truncation of the ChPT expansion commented 
above. Our estimates for the .fT-factors for p and a, given in the second part of Table HJ are 
an average of the results combining these different approaches, with errors taken to cover all 
the results. Let us note that the description of the a depends more strongly on the chiral 
loops, which are model- independent, and much less on the LECs. However, the dependence 
of the p resonance on the quark masses depends strongly on the values of the LECs. For that 
reason, for the central value of K q p we have only used the averaged results of the two best 
two-loop SU(2) ChPT fits in 25], which we consider to be the most reliable, in particular 
because they were fitted to three sets [53] of lattice calculations of the p mass dependence 
on the pion mass, which were consistent among themselves 4 , and because the resulting 
values of the LECs were more consistent with standard determinations and estimates. We 
25[ for details. For the strange quark mass dependence, we rely on the 



refer the reader to 



existing IAM one loop SU(3) calculations in Ref. [26], but including the uncertainties as just 
described above. 

It should be stressed that the quark mass dependences we find are significantly smaller 
than those given in Ref. jlij. In particular, in that reference a value of K s a = 0.54 is 
given compared to our —0.01 (c.f. Table |T]). The origin for this significant discre pan cy is 
the assumption about the quark structure of the a underlying the work of Ref. |l6j: the 
o was assumed to be an SU(3) singlet. In our case on the other hand the a emerges from 
nom-perturbative tttt interactions, which give only a very small dependence on the strange 
quark mass. 

As becomes clear from Table [U for all relevant quantities the variation with respect to the 
strange quark mass is smaller that the corresponding uncertainty. In addition, some quark 
mass variations driven by an external scale will lead to a relative change in the strange sector 
suppressed additionally by a factor m q /m s ~ 1/25. In what follows we will therefore only 
study the effect of a variation of the light quark masses on the NN potential. 



4 Let us nevertheless remark that there are other lattice calculations which are not quit e compatible with 
these three because their p masses fall systematically either above or below {55! them. The ones 
falling below are somewhat harder to accommodate within the IAM, as explained in |25| 
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III. IMPACT ON THE VARIATION OF THE TWO-NUCLEON POTENTIAL 



The changes in the meson masses cannot be directly connected to their impact on BBN. 
The quantity of relevance is the resulting variation of the two-nucleon (NN) interaction and, 
especially, its impact on nuclear binding energies. However, the connection of meson masses 



to the NN interaction is not clear cut. In Ref. 12] the phenomenological V18 interaction 
was used as a basis, where the insights of Refs. 17, [3] were translated into a variation of 



the model parameters. Although it provides some understanding of the sensitivity of the 
NN interaction on the variation of quark mass parameters, such a calculation is neither 



28 



an EFT approach is chosen; 



systematic nor complete. On the other hand, in Ref. 
however, there the input of the quark mass dependence of NN scattering lengths is taken 
from other sources. Our goal is to improve our understanding of the quark mass dependence 
of the NN observables using as input the .fT-factors given in Table [H 

As outlined in the introduction, we do not do the full dispersion theoretical treatment of 
Refs. 



35, 



in Ref. 



361 ]. but make connection to the NN force via the method of resonance saturation: 



31] explicit expressions are given that allow one to express the values of the NN 



contact terms in terms of meson masses and coupling constants. Thus, the quark mass 
dependences given above can be implemented straightforwardly. 



A. Quark mass dependence of the pion exchange contributions 

The long-range part of the NN potential up to next-to-next-to-leading order (N 2 LO) in 
the chiral expansion is driven by the exchange of one and two pions. In the exact isospin 
limit, the one-pion (I71") exchange potential at N 2 LO has the form 

V ^ ~ -4^2 \ gA ~ M ^ M *) Tl ' T2 g 2 + m 2 ' ^ ^ 

where cij denote the Pauli spin matrices, q = p ' — p is the nucleon momentum transfer and 
p (p ') refers to initial (final) nucleon momenta in the center-of-mass system. Further, F n 
and gA denote the pion decay and the nucleon axial coupling constants, respectively, while 
di$ is a low-energy constant from C n ^ that controls the leading contribution to Goldberger- 



Treiman discrepancy. Employing the spectral function regularization as detailed in Ref. 



59j, 
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the non-polynomial part of the two-pion (2tt) exchange potential has the form 



3^ 4 



647r 2 Fi 



£ A (<?) (^1 • 9^2 ■ q-ai-a 2 q 2 ) 



x 



A k {q) (a t ■ qa 2 ■ q - q 2 o x ■ a 2 ) , (3.2) 



with the loop functions L A (q) and A A (q) defined as 

72, ,2 i ^2„2 



L\q) = 6CA-2M n )^ln AuJ + q f + 2AqU * , uj = \/UAl + <f 2 , s = a/a 2 — 4M 2 
W V ; 2g 4M 2 (A 2 + g 2 ) V * V 



A A (g) = 9(k - 2M W ) — arctan g(A . (3.3) 

W 1 W 2g q 2 + 2AM 7r 1 ; 

Here and in what follows, the q are the LECs from the order-Q 2 (dimension two) pion- 



nucleon Lagrangian and A denotes the cutoff in the spectral representation, see Ref. [59 ]. 
In addition to the explicit M^-dependence, at N 2 LO one also needs to take into account the 
implicit one resulting from the chiral expansion of Qa and F n in the l7r-exchange potential 
in Eq. (13. ip . We use the NLO result for the chiral expansion of the pion decay constant as 
it is appropriate at the order we are working: 

F ^ f ( 1 + (i6^) l4+ °MJ ' (3 - 4) 

where F denotes the pion decay constant in the chiral limit. For the LEC I4 we adopt 
the value I4 = 4.3 from Ref. 60(. Using F n = 92.2 MeV, this leads to the chiral-limit value 
F = 86.2 MeV. Notice that the resulting NLO value for the if-factor K q FjT = 0.065 is slightly 
larger than the one given in Table HI 

Contrary to the pion decay constant and the nucleon mass, the chiral expansion for is 
known to converge rather slowly, see Fig. [TJ In particular, taking into account the leading 
(i.e. order 0(M 2 )) correction to its value at the chiral limit and adopting the value for 



the low-energy constant die — —1-76 GeV~ 2 obtained from the reaction ttN — > ttttN 62] 
leads to a very strong quark- mass dependence of g& near the physical point. On the other 
hand, lattice QCD calculations indicate that the behavior of gA with is rather flat. As 



discussed in Ref. 63J], such a flat behavior of g^, consistent with the lowest-mass lattice 
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0.5 



Q 



incomplete Q x ^ 
uncertainty in d J6 N 
variation of the lattice point \ 



J_ 



0,1 



0,2 

[GeV] 



0,3 



FIG. 1: Quark mass dependence of a a a t various orders in the chiral expansion. Also shown is the 



lowest mass lattice point from Ref. 



61[. The hatched band corresponds to the variation of die m 



the range di§ = —0.91 to —2.61 GeV , see Ref. 62|. The light shaded band results from a 10% 
variation of the lattice point used to fix the LEC of the operator in Eq. (|3.5p . 



data point from Ref. 6lj corresponding to M n = 353 MeV 5 can, in principle, be achieved 
at the two-loop level. In order to provide an accurate representation of the quark mass 
dependence of the l7r-exchange potential, we use in the present study the complete order-C} 3 
result accompanied by one order-Q 4 contact term, whose strength is adjusted to reproduce 
the lowest-mass lattice data from Ref. Q: 



9a 



9o 
+ 



f 4 T 1 f o , ox , M T 

1+ U* 6 -(^( 9 ° +(2+49 ° )ln K 



Mi 



(3 + 3g 2 - 4m c 3 + 8m c 4 ) M 3 + /3M 4 



(3.5) 



247rF 2 m 

Here, go, F and m refer to the chiral-limit values of the nucleon axial vector coupling, pion 
decay coup ling and the nucleon mass, respectively. We use the same values of the LECs q 
as in Ref. [63J namely c 3 = -4.7 GeV" 1 , c 4 = 3.5 GeV" 1 . 6 Further, M n = 138 MeV is 



We emphasize, however, that the convergence at such pion masses is problematic, as discussed in detail 



G4] 



in the review 

6 Notice that slightly different values of these LECs are adopted in the chiral NN potential. We have 
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the physical value of the pion mass in the isospin limit while (3 is a linear combination of 
the order-Q 4 LECs. We emphasize that the above expression does not correspond to the 
full order-Q 4 result since we do not include the order-Q 4 chiral logarithms. We have verified 
numerically that the effect of these logarithms is largely compensated by the /3-term when 
the later is tuned to reproduce the lattice point. Indeed, one observes that the solid line 
in Fig. [1] corresponding to the pion mass dependence of qa adopted in the present work is 



very similar to the more complete calculations of Ref. 



631 ] shown in Fig. 2 of that work. 



Further, Fig. [T] also shows the uncertainty associated with the variation of di§ in the range 



d 16 = -0.91 to -2.61 GeV" 2 



62] and the variation of the lattice point by 10%. 
The value of the nucleon mass in the chiral limit can be obtained from the order-Q 3 
expression 

m = m + 4 Cl M 2 + ^r 2 Ml + O (M 4 ) . (3.6) 

Using qa = 1.267 and c\ = —0.81 GeV -1 leads to m = 892 MeV. Note that the value of c\ 
used here is consistent with the small sigma-term, u^n = 45 MeV, cf. Eq. (I2.2p . We further 
emphasize that the resulting .fT-factor, K q N = 0.042, is consistent (within the error bars) 
with the value given in Eq. 



B. Quark mass dependence of the short-range terms 

The short-range potential in the 1 S'o and 3 Si — 3 D\ channels up to N 2 LO has the form 

^i5o° rt = C is , + C 1S o{p 2 + P i2 )i 

^35i rt = C 3S1 + C 3S1 (p 2 + P ' 2 ), 

V e f ort = C el p\ (3.7) 

where p = \p\, p' = 1^1 refer to the in-coming and out-going momenta in the center-of-mass 
system and Cj, Ci are M^-dependent coefficients 7 . The quark mass dependence of these 
operators can, in principle, be calculated straightforwardly in chiral perturbation theory 29 ]. 



verified that using these different values for C3.4 and re-adjusting the parameter j3 leads to a very similar 
Af^-dependence of gA- The induced difference in the two-nucleon observables is significantly beyond the 
theoretical uncertainty of our analysis. 
7 Of course, the LECs do not depend on the quark masses. The coefficients used here subsume the coefficients 
of the LO four-nucleon operators plus their first pion mass dependent corrections that are generated by 
the same operators times M%. 
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The problem is, however, that the coefficients in front of the contact operators oc M% are 



unknown. In 



291 ]. the corresponding LECs were estimated by means of naive dimensional 



analysis which, however, resulted in a very large theoretical uncertainty for two-nucleon 
observables. In order to avoid this difficulty, we follow here a completely different approach 

or contact operators. Instead, we make 



3l| and employ the M„--dependence for 



and refrain from doing an explicit chiral expansion 
use of resonance saturation of contact interactions 
the masses of the heavy mesons discussed in Sec. Ill Bl 

Resonance saturation for contact NN operators up to N 2 LO is discussed in detail in 



Ref. [3JJ. In that work strongly reduced values of the LECs q were adopted in order 

(2tt) 



59|, 



to circumvent a very strong attraction resulting from the isoscalar central two-pion 
exchange potential calculated using dimensional regularization. As discussed in Ref. 
the strong attraction in the 27r-exchange potential at intermediate and shorter distances can 
be traced back to the large-mass components in the spectrum which cannot be described 
reliably within the framework of chiral EFT. In the chiral potentials of Refs. 65|, |66(, the 
unphysical high-mass components in the two-pion exchange spectrum are cut off by means 
of the spectral function regularization. We now repeat the analysis of Ref. 31] for the actual 
version of the chiral N 2 LO potentials and test the validity and accuracy of the resonance 
saturation hypothesis. 



Here and in what follows, we consider the Bonn-B [67| potential as a typical represen- 
tative of one-boson-exchange (OBE) models. Its long-range part is given by l7r-exchange 
whereas shorter-distance physics is expressed as a sum over contributions from the exchange 
of heavier mesons. For nucleon momentum transfers well below the masses of the exchanged 
mesons, one can interpret such exchange diagrams in terms of local contact operators with an 
increasing number of derivatives (momentum insertions). The LECs accompanying the re- 
sulting contact interactions are then expressed in terms of the meson masses, meson-nucleon 
coupling constants and the corresponding form-factors. In order to allow for a meaningful 
comparison between the chiral and OBE potentials, one needs to properly account for the 
chiral 27r-exchange potential, which contributes to the chiral potential but is absent in the 



OBE models. Here we follow the strategy of Ref. 



3JJ and expand the 27r-exchange potential 



of Eq. ( 13. 2 p in powers of momenta. This induced contributions to the LECs entering the 
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1 5'o and 3 Si — 3 D\ partial waves read: 



ISO 



00 150 



-8g\ + 4gj + 1) M 2 X /A 2 - 4M 2 
247rFiA 



A2 _ 4M 2 ( ( _ 8% 4 + 17g 2 A + 2 ) A 2 + 2 (8g% - Ag\ - 1) M 2 

144vrFiA 3 



(8g\ - - 1) M^ Jh* - AMI 
87rFiA 



v A 2 - AMI ((40^1 - 17s 2 - 2) A 2 + (-16^ + 8g\ + 2) M 2 
351 ~ 48vrF4A 3 ' 

60^ = -^-= ^JL, (3.8) 

£l 4V27TF4A V ; 

at NLO and 

y(3) 3^M 3 (2c 1 -c 3 )(2M 7r -A) 



^ 15 ° " 4FM 

7T 

S 2 M 7r (2M 7r - A) ( A 2 (-10 Cl + llc 3 - 4c 4 ) + 4M 2 (2 Cl - c 3 ) + 2AM 7r (2c 1 - c 3 ] 

SC[% = * = 

150 16F T 4 A 3 

sd (3) _ 3s 2 M 3 (2c 1 -c 3 )(2M,-A) 

7T 

C3, 



S 2 M,(2M, - A) A 2 (-10 Cl + llc 3 - 4c 4 ) + 4M 2 (2 Cl - c 3 ) + 2AM,(2 Cl 
^gi = * 



16FiA 3 



y(3) c 4 g\M n (A - 2M W ) 



5C w = ^ — (3.9) 

at N 2 LO. In the limit A — > oo corresponding to dimensional regularization, the above ex- 



pressions agree with the ones given in Ref. 



311 ] . The size of the 27r-exchange-induced con- 



tributions to the LECs for the two extreme values of the spectral function cutoff A can be 
read off Table [TTJ 

After these preparations, we are now in the position to test the resonance saturation 
hypothesis for Ci and Cj. The contributions of the various mesons to the contact oper- 
ators can be obtained by carrying out partial wave decomposition of the expressions for 



the boson exchange contributions given in Ref. 31] and expanding the results in powers of 



momenta. In Tables II III and IIVt the 27r-exchange-corrected values of these LECs for five 
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LEC 


Q 2 , A = 500 MeV 


Q 2 , A = 700 MeV 


Q 3 , A = 500 MeV 


Q 3 , A = 700 MeV 


5C\so 


-0.004 


-0.005 


-0.004 


-0.005 




-0.534 


-0.592 


-0.365 


-0.500 


SC 3S i 


0.013 


0.014 


-0.004 


-0.005 




0.594 


0.663 


-0.365 


-0.500 




-0.178 


-0.196 


0.170 


0.229 



TABLE II: Contributions to the LECs Q and Cj induced by the NLO and N 2 LO 27r-exchange 
potential. The Ci are in 10 4 GeV" 2 and the d in 10 4 GeV~ 4 . 





fit 1 fit 2 fit 3 fit 4 fit 5 


Bonn-B 


Ciso + ^llo 
Ciso + <5Ciso 

C3SI + ^li 
Cssi + ^ 
C 6l + Scf? 


-0.161 -0.066 -0.095 -0.161 -0.111 
0.974 1.574 1.457 1.008 1.386 

-0.169 -0.136 -0.135 -0.167 -0.135 
0.356 0.256 0.231 0.280 0.221 

-0.390 -0.332 -0.325 -0.373 -0.321 


-0.117 
1.276 

-0.101 
0.660 

-0.410 



TABLE III: LECs C% and Ci from the NLO chiral potential for different cutoff combinations (fits 1 
to 5 as defined in Eq. (|3.10p ) corrected by the induced contributions of the 2-7r-exchange potential. 
Also shown are values resulting from resonance saturation using the Bonn-B model (last column). 
The Ci are in 10 4 GeV" 2 and the d in 10 4 GeV" 4 . 



cutoff combinations 

NLO : {A, A} = {400, 500}, {550, 500}, {550, 600}, {400, 700}, {550, 700} , 
N 2 LO : {A, A} = {450, 500}, {600, 500}, {550, 600}, {450, 700}, {600, 700} , 

(3.10) 

are listed together with the values resulting from resonance saturation based on the Bonn-B 
model. For the considered cutoff combinations, all LECs are reasonably well described in 
terms of resonance saturation. The agreement is better at N 2 LO, most notably for the LEC 
C^si- The estimations based on resonance saturation yield the numbers which are typically 
within 20 - 30% of the true values except for the LO LECs Q which appear to be somewhat 
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fit 1 fit 2 fit 3 fit 4 fit 5 


Bonn-B 


Ciso + SC[% + 5C[% 

Ciso + sd 2 i + scfi 

loU loll 

C3SI + + ^^Sl 

<? 3 si + ScQx + Sci% 
C £ i + 6C® + <5C e ( ? 


-0.161 -0.116 -0.159 -0.163 -0.161 
1.164 1.368 1.243 1.321 1.321 

-0.162 -0.127 -0.137 -0.164 -0.130 
0.574 0.423 0.523 0.720 0.568 

-0.425 -0.363 -0.395 -0.467 -0.409 


-0.117 
1.276 

-0.101 
0.660 

-0.410 



TABLE IV: LECs Ci and Cj from the N 2 LO chiral potential for different cutoff combinations (fits 1 
to 5 as defined in Eq. (I3.10p ) corrected by the induced contributions of the 27r-exchange potential. 
Also shown are values resulting from resonance saturation using the Bonn-B model (last column). 
The Ci are in 10 4 GeV" 2 and the C l in 10 4 GeV" 4 . 



more strongly underestimated. For the sake of completeness, we also list in Table IV 
contributions from individual mesons exchanges in the Bonn-B model, see also Ref. [31 1 



the 



LEC 


V 





5 


ijj 


P 


sum 


/^res 
ISO 


0.000 


— 0.392 


-0.023 


0.287 


0.011 


-0.117 


fires 
°1S0 


0.033 


1.513 


0.036 


-0.560 


0.254 


1.276 


°351 


0.000 


-0.424 


0.070 


0.287 


-0.034 


-0.101 


fires 
°3S1 


-0.011 


1.030 


-0.108 


-0.777 


0.526 


0.660 


/~tres 


-0.032 


0.000 


0.000 


0.077 


-0.455 


-0.410 



TABLE V: Contributions of the various boson exchanges to the LECs for the Bonn-B potential 
and the corresponding sum. The Ci are in 10 4 GeV -2 and the Ci in 10 4 GeV -4 . 



The observed reasonably good representation of the LECs accompanying the short-range 
operators in terms of heavy-meson exchanges justifies modelling the chiral extrapolations for 
the corresponding operators in terms of quark/pion mass dependence of the heavy mesons as 
discussed in Section UTBl More precisely, we assume that the resonance saturation hypothesis 
remains valid at unphysical values of the quark/pion masses, that is 

X Z (M V ) + SXriM,) = XJ(M n ) + X%(M n ) + Xf(M v ) + X}^ . (3.11) 
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FIG. 2: Various contributions to the quark mass dependence of the LECs C351 (left panel), C^si 
(middle panel) and C e \ (right panel) for fit 1 at N 2 LO as discussed in the text. Here, m g o denotes 
the physical value of the light quark mass. 



Here X stays for C, C and I = {ISO, 3SI, el} and 5X I (M n ) = SX\ 2) (M n ) at NLO and 
6Xi(M v ) = 8Xf\M % ) + 8Xf ) {M„) at N 2 LO. For the resonance contributions, we take into 
account the quark mass dependence of the a- and p-meson masses as given in Table IT] and 

nn 

assume K% = K q for the a;- meson [13, HSJ. Neglecting the quark mass dependence of the 
77- and 5-mesons is justified by their small contributions to the LECs. The resulting error is 
expected to be well below the theoretical uncertainty of our analysis. Notice that in X° ,p,u} , 
we also take into account the quark mass dependence of the nucleon mass. The last term on 
the right-hand side of Eq. ( 13. lip . XJ est , denotes the contributions to the LECs not related 
to the heavy-boson exchanges. The (unknown) M^-dependence of Xj est is neglected in the 
present work. This, of course, only makes sense if Xj cst is small compared to the Xj, i.e. if 
the LECs are well described in terms of resonance saturation. This is indeed the case for 
both the NLO and N 2 LO potentials. As a representative example, we show in Fig. [5] the 
individual contributions to the quark mass dependence of 6*351, C351 and C £l resulting from 
Eq. (13. lip for fit 1 at N 2 LO. While strong cancellations between the a- and ^-contributions 
are observed for C^si and C351, the LEC C& is largely saturated by the p- meson. Notice 
that as a result of the cancellations, there is a sizeable uncertainty in the m g -dependence of 
C351 associated with the non-resonance contribution of the last term in Eq. ( 13. lip . 
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C. Quark mass dependence of the S-wave iViV observables 



Having specified the quark mass dependence of the NN potential, we now turn to the 
chiral extrapolations of two-nucleon S-wave observables. We calculate the NN phase shifts 
and mixing angles by solving the nonrelativistic 8 Lippmann-Schwinger equation in the partial 
wave basis 

w,p) = w,p) %$- v iwy) p *_^ +i / ti(p",p) , (3-i2) 

as rj — > + . The relation between the on-shell S- and T-matrices is given by 

s 'n (p) = s i'i - ^P m T n (p) • ( 3 - 13 ) 

Thus, the quark-/pion-mass dependence in the observables emerges from both the nucleon 
mass entering the NN Green's function and the potential. 

Our results for the chiral extrapolation of the deuteron binding energy and the inverse 
S-wave scattering lengths at NLO and N 2 LO are visualised in Figs. [3] and HI In these figures, 
m q o denotes the physical value of the light quark mass. In these calculations, we also took 
into account the implicit quark mass dependence in the two-pion exchange potential induced 
by Qa an d F n which is, strictly speaking, a higher-order effect. We will comment on the size 
of these higher-order contributions in the next section. Also, as already explained above, we 
follow here our general strategy and use the most accurate available information regarding 
the m q -dependence of m and, especially, of gA coming, in particular, from lattice QCD 
simulations rather than sticking to the strict chiral expansion at a given order. 

We observe the opposite trends in the l So and 3 S\ channels when changing the value of 
the quark mass. In particular, the interaction between the nucleons is found to become more 
attractive in the 1 Sq channel with increasing the light quark masses while more repulsive 
in the 3 S*i partial wave. The obtained results do not exclude the possibility of a bound 
spin-singlet state at sufficiently large quark masses. The deuteron is found to remain bound 
for all values of the quark masses considered. Notice further that our results indicate that 

n n 

the infrared limit cycle proposed in Ref. [68] (see also Ref. [69[) is unlikely to emerge in the 



Relativistic corrections to the two-nucleon Green's function need to be taken into account starting from 
N 3 LO which is beyond the scope of the present work. 
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FIG. 3: Quark mass dependence of the deuteron binding energy at NLO (light-shaded band) and 
N 2 LO (dark-shaded band). The bands correspond to the cutoff variation as discussed in the text. 
The solid dot shows the deuteron binding energy at the physical quark mass. 



range of the quark masses considered in the present analysis. A detailed comparison of our 
findings with the available calculations will be presented in the next section. 

Let us now address the theoretical uncertainty of our calculations. It is comforting to 
see that the results for the quark mass dependence of the deuteron binding energy and the 
S-wave scattering lengths calculated at NLO and N 2 LO are consistent with each other. The 
NLO and N 2 LO bands resulting from the cutoff variation as described above overlap except 
for large quark masses in the spin-triplet channel. This is in line with the observation 
that the cutoff variation at NLO underestimates the true theoretical uncertainty at this 
order since the width of the bands at both NLO and N 2 LO measures the impact of the 
neglected order-Q 4 contact interactions. A more complete discussion of the theoretical 
uncertainties of the calculated NN observables will be given in Section IIVI Notice further 
that the chiral extrapolations become rather uncertain for small quark masses which, at first 
sight, might appear counterintuitive. This, however, has to be expected given that the LECs 
accompanying contact interactions are fit to experimental data which obviously correspond 
to the physical quark masses. In addition, one should also keep in mind that we do not 



rely here on the chiral expansion of the short-range forces, contrary to Refs. 



29. 



34|. Our 
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FIG. 4: Quark mass dependence of the inverse S-wave scattering lengths at NLO (light-shaded 
band) and N 2 LO (dark-shaded band). The bands correspond to the cutoff variation as discussed 
in the text. The solid dots show the inverse scattering lengths at the physical quark mass. 



approach, utilizing resonance saturation and the i^-factors for heavy-meson masses, cannot 
be expected to yield reliable results at low quark masses where short-range contributions 
non-analytic in quark masses, which are not explicitly taken into account in our calculations, 
should play an important role. 

Finally, we list in Table |VT] the individual contributions of various mechanisms to the 
dimensionless quantities K^ ent and K\ x with / = {150, 351} defined as 



K. 



SE deut /5m q 



deut 



K q 



5ai/5m q 



(3.14) 



E deut /m q ' °> ai/m q 

We observe a reasonable convergence pattern in the triplet channel with the main effect 
coming from the LO terms in the potential and the contributions due to NLO+NNLO terms 
and the quark mass dependence of the nucleon mass being considerably smaller. The much 
larger uncertainty in the singlet channel can be explained by the known feature that the one- 
pion exchange yields only a minor contribution to the 1 Sq phase shift. The dominant effects 
emerge from two-pion exchange and shorter-range terms whose quark mass dependence is 
less constrained than the one associated with the longest-range one-pion exchange potential. 
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K q 

a 


ISO 


if 9 

o 


3S1 


K q 

A deut 






NLO 


N 2 LO 


NLO 


N 2 LO 


NLO 


N 2 LO 


V 1tt + Kont 

+ m (LS eq.) 


2 +o.i 

i,J -0.3 


-0-8^ 

9 q+1.6 

z,a -1.5 


n 36 +0 - 09 

24+0.O6 

U -^-0.34 

n i o +0.05 
U - i<5 -0.33 


54 +0 ' 00 

U - J -0.05 

43+o.os 

-0.11 

n or, +0.08 

u - JZ -o.io 


n o 7 +o.06 

u,uu -0.13 

_ 41 +0.76 _ 
-0.13 


-, 9 n+0.12 

1.28 _ Q 

1 1 1 + - 27 
1,11 -0.19 

0.8633 



TABLE VI: Various contributions to K q a igo , if^ 3S1 and ^ eut - The numbers correspond to the 
third cutoff combination in Eq. (|3.1U|) whiie the errors resuit from the cutoff variations. 



IV. DISCUSSION 



We are now in the position to discuss in some detail the theoretical uncertainty of our 
calculations. Its main sources are due to 

1. The uncertainty associated with the chiral extrapolation the nucleon mass m as well as 
the axial coupling g& and the pion decay constant F n , which impact the m^-dependence 
of the long-range interactions. 

2. The uncertainty due to truncating the chiral expansion of the potential at N 2 LO. 

3. The uncertainty associated with the resonance saturation hypothesis for short-range 
operators and the employed quark mass dependence of the heavy-meson properties. 

The chiral dependence of is well reproduced at the leading-loop order in ChPT, so 
that the associated uncertainty has a much smaller impact on the two-nucleon observables 
considered in this work as compared to other sources. On the other hand, the chiral expan- 
sion of qa is known to converge slowly. As explained in Section IIII A[ we incorporate the 
order-Q 4 counterterm and use the lattice QCD result at M n = 353 MeV as an input to tune 
the corresponding LEC. This allows us to obtain a realistic description of the quark mass 
dependence of g&. We check the robustness of this procedure by allowing for a 10% variation 
of the lattice point. As visualized in the left panel of Fig. El this induces a shift in the binding 
energies which is considerably smaller than our estimated theoretical uncertainty, cf. Fig. |3J 
On the other hand, the uncertainty in the determination of the LEC d\%, die — —0.91 to 



-2.61 GeV- 2 



621 ]. leads to a sizeable spread which is comparable with the one emerging 
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g A : uncertainty in d [6 

g A : 10% variation of the lattice point 

,3 




> 

w 



uncertainty in K 
uncertainty in K 
uncertainty in K 
NNLO, fit 3 



m / m . 

q qO 




FIG. 5: Quark mass dependence of the deuteron binding energy at NNLO, fit 3 (solid lines). 
Left panel: theoretical uncertainty associated with the quark mass dependence of the long-range 
interactions and the nucleon mass. The hatched band corresponds to the variation of d\Q in 



the range d 16 = -0.91 to -2.61 GeV" 2 (see 



62]]). The light shaded band results from a 10% 



variation of the lattice point used to fix the order-Q 4 counter term in Eq. (|3.5p . Finally, the long- 
dashed-dotted line shows the effect of neglecting the quark mass dependence of gA and F n in the 
27r-exchange potential while the short-dashed-dotted line shows the effect of using the order-Q 2 
rather than order-Q 3 expression (]3.6|) for the chiral extrapolation of the nucleon mass. Right panel: 
theoretical uncertainty induced by the errors in the quark mass dependence of the heavy meson 
masses according to Table HI 



from the cutoff variation. This can be expected since the value of die influences the shape of 
the quark mass dependence of (larger in magnitude values lead to a more flat behavior), 
see Fig. [TJ 

The uncertainty due to truncating the chiral expansion for the potential at N 2 LO was 
already roughly estimated by the cutoff variations, see Fig. [3j As an additional check, we 
calculated the impact of the M^-dependence of g& and F n in the 27r-exchange potential 
which nominally starts to contribute at N 3 LO (i.e. order Q A ) which is beyond the accuracy 
of this work. The size of this effect is given by the difference between the solid and long- 
dashed-dotted lines in the left panel of Fig. and is indeed within the estimated theoretical 
accuracy at N 2 LO. Similarly, it is, strictly speaking, sufficient to use the order-Q 2 rather 
than order-Q expression for the M w - dependence of the nucleon mass at N 2 LO. The induced 
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difference agrees well with the estimations based on dimensional arguments and is within 
the accuracy of our calculation, see the short-dashed-dotted line in the left panel Fig. [5j 

Presumably, the most important source of uncertainty is due to the quark mass depen- 
dence of the contact interactions. While resonance saturation itself seems, at least in prin- 
ciple, to provide a fairly accurate way to relate the chiral extrapolations of the short-range 
terms to the properties of heavy mesons which are easier accessible theoretically, it is diffi- 
cult to estimate the theoretical uncertainty associated with this procedure. We, therefore, 
restrict ourselves to propagating the errors in the i^-factors for the heavy meson masses, 
see Table [J through our analysis. These errors turn out to be strongly magnified due to 
the large cancellations between the contributions of the p and u mesons, see Fig. [2J The 
resulting uncertainty in the deuteron binding energy appears to be comparable to the one 
emerging from the cutoff variation and the chiral extrapolations of Qa, see Fig. |5j We further 
emphasize that using the linear approximation in terms of the i^-factors for the quark mass 
dependence of the heavy mesons is yet another approximation (if one goes sufficiently far 
away from the physical point). It is, however, irrelevant as long as one is only interested in 
the .fT-factors and can be easily avoided if necessary. 

The final results for the i^-factors of the deuteron binding energy and the corresponding 
scattering lengths read 



K, iso = 2 3 _ 1 g > K, 3Si = °- 32 _ i 8 ' ^cut = _0 - 86 _o 50 ' 



(4.1) 



where the numbers are given for the third cutoff combination at N 2 LO and the central values 
of K% listed in Table HI The theoretical uncertainties due to truncating higher order terms 
(estimated by the cutoff variation), the uncertainty in d 16 , the lattice calculation of g^, and 
the errors in K% are added in quadrature. 

The results given in Eq. (14. ip can be compared with the ones of Ref. j^, carried out at 
NLO in chiral EFT 

^a,iso = 5±5, ^38! = 1.1 ±0.6, Kl ut = -2.8±1.2, (4.2) 

where the numbers are inferred from Figs. 11 and 12 of this work. A more conservative 
error estimation taking into account a larger variation in the LEC die an d in the quark mass 
dependence of the lowest-order spin-triplet contact interaction was performed in Ref. [70]] 
leading to 

K 3Si = LI ± 0.9 , Kl nt = -2.9 ± 1.8 . (4.3) 
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As already 



one of 
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pointed out before, the major differences between the present analysis and the 



70| are in using a more realistic result for the chiral expansion of g^, employing 



resonance saturation to describe the quark mass dependence of contact interactions and 
extending the calculations to N 2 LO. These improvements have allowed us to determine the 
values for the .fT-factors with higher accuracy. 

The results for both 1 5'n and 3 Si channels given above are consistent with the chiral EFT 
values calculated in Ref. 



K 1S0 = 2-4 ± 3.0 



7l| using the Kaplan-Savage- Wise approach, 

K 3Si = 3.0 ±3.5, i^ eut = -7±6 



(4.4) 



where the numbers correspond to Figs. 1, 2 and 4 (rj — 1/5) of that work. 

More recently, attempts have been made to combine chiral EFT with lattice-QCD cal- 
culations. In particular, the NPLQCD collaboration has determined the regions for the 
S-wave scattering lengths consistent with their lattice data a\so = 0.63 ± 0.50(5-10) fm and 
a 3S1 = 0.63 ± 0.74(5-9) fm obtained at M n = 353.7 ± 2.1 MeV |z2|. Assuming the validity 
of the employed chiral EFT frameworks in the considered range of pion masses, they deter- 
mined the following constraints for the allowed regions of the f^-factors (the numbers below 
are extracted from Figs. 3, 4 of Ref. |72j ) : 



K q < — 4 1 1 K q > ? 

-"■a, ISO ~ * U ^a, ISO ~ Z 



based on the Weinberg approach and 



K 1S0 > 6 



-5 <K, 3Si^-0-2 U 0.3<i^3 S1 <9, 



(4.5) 



(4.6) 



using the chiral EFT formulation of Ref. 



73|. For a more recent extrapolation of the 

n 

NPLQCD numbers to physical quark masses see Ref. [74J- Very recently, a similar analysis 



has been carried out using the framework of chiral EFT with dibaryon fields 



33| yielding 



K q 3S1 ~ —0.4, see Fig. 9 of that work, and a positive value for K q a 1S0 (which we were 
unable to infer from that figure). While this result for K q 3S1 disagrees with our findings, it 



is difficult to draw conclusions since Ref. 



331 ] does not provide an estimate of the theoretical 



uncertainty associated with using the lattice-QCD results at large values of the pion mass 



72J given 



as input in the calculations. The same comment also applies to the results of Ref. 
in Eqs. fT43|) . (T4l| as well as the results of Ref. Q. 

Last but not least, it is comforting to see that our results are in a good agreement with 



the ones of Ref. Il2j, where the value i^ eut = —0.75, —0.84 and —1.39 is reported for three 
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different models of the two-nucleon potentials. Even more important is, however, that we 
are able to carefully estimate the theoretical uncertainty for this quantity. 



V. CONSEQUENCES FOR HEAVIER NUCLEI 



So far we focused on the two-nucleon system, however, the quark mass dependences of 
3 He and 4 He are also relevant for BBN. In order to estimate the impact of the quark mass 
dependences of the deuteron binding energy and the 1 Sq scattering length on BBN, we here 



use the methods of Ref. 



28j - actually, what we have provided in the previous sections is 



an update of the input used in Ref. |28|, which basically came from Ref. [29]. 

The quark mass dependences of the binding energies of 3 He and 4 He can be calculated 
from 



a, ISO 
a, lSO-^He 



K q K 



+ A deut A ^ 



dcut 
4 He > 



(5.1) 



where the K% Re = xj AA He (<5£U He / 5x) were defined in analogy to the quantities given above. 
In Ref. [28( the coefficients were calculated from pionless EFT. They read, including the 



uncertainties quoted in Ref. 



28) 



= 0.12 ±0.01 , 
A 4 a ^f° = 0.037 ±0.011 , 



Ki^ = 1.41 ± 0.01 ; 
Xfeut = g 74 ± Q.22 ; 



(5.2) 



From this we get using the uncertainties for the A-factors as given in Eq. (14. ip 



K$ Re = -0.94 ± 0.75 



K! Re = -0.55 ± 0.42 , 



(5.3) 



where the uncertainties of Ref. 



28 



and those quoted in Eq. (14.1 1) were added in quadrature. 
Note that there has been a recent lattice study of the deuteron, 3 He and 4 He at a pion mass 



of 510 MeV and various lattice sizes 



75| . As this pion mass is rather large, we refrain from 



trying to extract the corresponding A-factors from that work. In the future, when such 
studies become available at lower values of M n , they will provide valuable constraints on the 
quark mass dependence of nuclear binding energies. Note further that the A-factor for 4 He 
is consistent for the central value obtained from a recent nuclear lattice calculation using 



the same scattering lengths A-factors, Kl He = —0.32 [7] 
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X 


d 


3 He 


4 He 


6 Li 


7 Li 


a s 


-0.39 


0.17 


0.01 


-0.38 


2.64 


-E>deut 


-2.91 


-2.08 


0.67 


-6.57 


9.44 


-Btrit 


-0.27 


-2.36 


0.01 


-0.26 


-3.84 




-2.38 


3.85 


0.01 


-5.72 


-8.27 


#4 He 


-0.03 


-0.84 


0.00 


-69.8 


-57.4 


Beu 


0.00 


0.00 


0.00 


78.9 


0.00 




0.03 


0.01 


0.00 


0.02 


-25.1 


Br Be 


0.00 


0.00 


0.00 


0.00 


99.1 


T 


0.41 


0.14 


0.72 


1.36 


0.43 



TABLE VII: BBN response matrix d In Y a /d In X at baryon-to-photon ratio r\ = 6.19 x 10~ 10 . The 
Y a are the number ratios of primordial isotope abundances to hydrogen, except for 4 He which is 
the mass ratio 4 He/H. 

VI. IMPACT ON BBN 

In Table I VII I we present our calculated BBN response matrix. The quantities presented 

are the linear dependences of calculated primordial abundances Y a to small variations of 

nuclear binding energies and scattering lengths X, d In Y a /d In X. They were obtained using 

I I 

the methods and codes (modified from the publicly available Kawano code [761 ]) described 
in Updated reaction rates are taken from Refs. 77-8oj], see [kJ for details. Unlike in 



previous studies, we have separated the effect of -Bdcut from the singlet scattering length a 



The scattering length affects the rate of the reaction n(p, d)j via the equation (see, e.g. [8l|) 

(av) ~ (5 dcut ) 5/2 A, 

where e v is the position of the singlet virtual level. The baryon-to-photon ratio rj = 
6.19 (15) x lCr 10 is taken from the latest WMAP7 data 3]- 

Final sensitivities of primordial abundances Y a to variation of m q , are obtained by com- 
bining the results in Table IVHI with Eqs. (14.11 15. 3ft using 

SlnY a _^ dlnY a 

6\nm q ^ dlxiXi Xr K ' J 
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Ya 


obs. 


calc. 


5lnY a /Sln m q 


Sm q /m q 


[deut/H] 


2.82(21) x 1(T 5 


2.49(17) x 1(T 5 


3.9(3.4) 


0.034 (42) 


4 He (Y p ) 


0.249 (9) 


0.2486(2) 


-0.56(34) 


-0.003 (65) 



TABLE VIII: Extracted values of quark mass variation, 5m q /m q , during BBN from comparison of 
observed and calculated primordial abundances, Y a . Observed values are taken from the Particle 



Data Group review [83] and calculated values from Ref. [841 ] . 

Error estimates must be performed carefully, since the K q x are derived from the same sources 
via Eq. ( 15. ip . We have taken the correlations in errors into account when deriving our final 
sensitivities. The uncertainties of Eq. ( 15. 2 p are also correlated to some extent, but since 
they are small anyway we may neglect this correlation. The final sensitivities of primordial 
abundances to quark mass variation are given in Table IVIIIl along with final values of quark 
mass variation at the time of big bang nucleosynthesis extracted from comparison of observed 
and calculated primordial abundances. 

We see from Table IVTlTl that the deuterium and 4 He data give consistent limits on 5m q /m q . 
Taking a weighted average of the two results gives Sm q /m q = 0.02 ± 0.04, our final result. 

VII. EFFECT OF THE NEUTRON LIFETIME 

Our limit 5m q /m a = 0.02 ± 0.04 at first appears much weaker than the limit derived 
by Bedaque et al. [28(, who obtained —1% < 5m q /m q < 0.7%, although our input for the 
two-nucleon parameters is more accurate. The origin of the difference is that, contrary to 
this work as well as previous works [lO, 12, 13], in Ref. j^sj a quark mass dependence of the 



neutron lifetime, r, was included. Since essentially all free neutrons at the onset of BBN 
end up in 4 He nuclei, and changing r changes the neutron-to-proton number ratio at BBN, 
the 4 He abundance is quite sensitive to the neutron decay rate. Therefore it is worthwhile 
to examine the assumptions made in Ref. [28j in more detail. 

Neutron beta decay becomes possible as a consequence of a non-vanishing proton-neutron 
mass difference, which is non-zero due to an apparent violation of the isospin symmetry in 
the Standard Model caused by m u ^ rrid and electromagnetic effects driven by different 
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quark charges. One finds 

Am N = m n -m p = Amf + Am c N m = 2 MeV - 0.7 MeV, (7.1) 



with an uncertainty of 0.3 MeV for the individual contributions 85J 9 . Thus, in order to 
quantify how Amjv and thus r changes with the quark masses, an assumption has to be 
made on how the other Standard Model parameters change — as we will see of particular 
importance is the change of the electron mass relative to the quark masses. 
The neutron width F ~ 1/r can be written as 



(G F cos6 c ) 2 5 . 2 \ / Am 



_^ mt(1+3A)f[ =^ ) (7 . 2) 

where Gp is the Fermi constant, 6c is the Cabibbo angle, m e is the electron mass, and g& is 
the nucleon axial decay constant. The function /(AfiWm e ) describes the phase space and 



Coulomb attraction; an explicit form is presented in 



28] 



In order to proceed, Ref. [28j| made the assumption that when m q changes, m u /m^ as 
well as all other Standard Model parameters stay constant. This clearly introduces some 
model-dependence. Based on this assumption one gets a very strong sensitivity of V to m q 
via the factor /: 



SlnT Slnf(Am N /m e ) f 



S In m q S In m g f 



p(AWm e ) \ 4 
\ o In m q J 



where /'// is numerically determined to be 2.57 at the physical value of Am^jm e . Thus, 
the large sensitivity to the variation of the quark mass, found in Ref. comes directly 
from mo del- dependent assumptions, which allow one to write 

5{ArriM /m e ) / 5 lnm g = Am% T /m e = 4 . 

One the other hand, had one assumed that m u — is independent of the quark mass, the 
value found for b\\iY j b\wm q would be significantly smaller. 

In general, how possible changes of fundamental parameters are interrelated depends 
on the model assumed for the physics beyond the Standard Model. In fact, the relation 
mu/rnd = constant emerges naturally from a scenario where all elementary particle masses 
are proportional to a varying Higgs vacuum expectation value, v (relative to Aq CD ), but 



9 The more recent extraction of Ref. [86| finds Am^J 11 = — 1.5 ± 0.5 MeV — consistent within uncertainties. 
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_ < 


v m e J 



3.6 ±1.5. (7.4) 



the gauge and Yukawa couplings are independent of it. However, in this case m e and the 

mass of the weak gauge boson, mw, are also proportional to v. One finds, therefore, that 

Am^ r /m e has no dependence on v, and the sensitivity of / comes from the variation in m e 

relative to the electromagnetic component of Am^: 

S\nf(Am N /m e ) _ f_ 
Slnv f 

Under our assumption that Slnrriw = S\nm e = 5lnm q — Shxv (i.e. all masses vary with 
the Higgs VEV), and noting that Gf ~ l/ m iy; from (17. 2p we obtain 

^ = 1 + : ^ 2 5 !y l + 6 ±I^M = 4,8 ± 1.5 , (7.5) 
6 In v 1 + 6g A omv omv 

where we used the quark mass dependence of g\ discussed in Sec. IIII Al which gives 

8\ng 2 A 



5lnv 



0.2 ±0.1. (7.6) 



The final sensitivity to neutron decay rate, 8\nr / 5\nv = —4.8, is reduced by around a 
factor of two compared with j^j, but it is still very large even within our model. Note, 



Shir/ Shxv = —4.9, may also be obtained using Table IV of Dent et al. 10] under the same 
assumptions, and a value of —4.8 can be extracted from Eq. (7) of Ref. [87J. Multiplying 
Slnr/Slnv by the sensitivity coefficients 5lnY a /5lnT presented in Table rvTTI and adding 
the binding energy and scattering length sensitivities presented in Table IVIIII (using our 
assumption 5h\m q = 5\nv), we obtain 

= 1.9 (3.4) , = -4.0 (0.3) . (7.7) 

o In v o In v 

In this model, the r sensitivity reduces the deuterium sensitivity to v by a factor of two, 
and entirely dominates the 4 He sensitivity. The final limits on Higgs VEV variation from 
deuterium and helium abundances are dv/v = 0.07±0.13 and 5v/v = 0.000±0.009, respec- 
tively. However, we should put here the disclaimer that, although we included the effects 
of isospin violation in the evaluation of the neutron beta decay, all few-nucleon calculations 
were done imposing isospin symmetry. One expects these effects to be small compared to 
that from neutron beta decay; still, for consistency such effects will have to be included in 
the future. 

Therefore, imposing the model that all masses scale with v and including the variation 
of r, then as expected the results are entirely dominated by the 4 He data and the lim- 
its are rather tight, |5t>/t>| < 0.9%. Finally, we note that this result does not include all 
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possible mass variations: for example, in Ref. [10| the dependences of primordial abun- 
dances on m e are calculated assuming r constant (that is, the m e dependence that does not 
come via r variation). In comparison with the dependence from r these are rather small, 
5lnY a /5lnm e = —0.16 and —0.71 for deuterium and 4 He, respectively. Including this effect 
would marginally tighten our limits on Sv/v. 

VIII. SUMMARY 

We have presented a systematic study of the impact of quark-mass variations on prop- 
erties of the two-nucleon system based on chiral perturbation theory combined with non- 
perturbative techniques. Since the approach is based on an effective field theory a reliable 
error estimate becomes feasible — a clear advantage compared to purely phenomenological 



studies as reported, e.g., in Refs. [16H21|. We include the uncertainties from the quark mass 
dependence of meson masses as well as those from the treatment of the NN interaction de- 
rived from chiral perturabtion theory. Here the heavy (heavier than the pion) meson masses 
enter in the expressions for four-nucleon contact terms with strengths fixed from a resonance 
saturation hypothesis — the uncertainty induced by this method is not included. However, 
support for this procedure comes from the apparent quantitative success of the resonance 
saturation for four-nucleon contact interactions. Especially we find 

K, 1B0 = 2-3 + ^ , Kl 3S1 = 0.32+°^ , Kl nt = -0.86^ Q f Q , (8.1) 

where the uncertainties are significantly reduced compared to the numbers derived from 



earlier studies within effective field theory 
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30] . The numbers for the A"-factors presented 



in this work are the necessary input for studies that address the quark mass dependence of 
nuclear properties. These studies allow one to quantify how much fine-tuning is necessary 
amongst the Standard Model parameters to allow life to develop — for the most recent 
developments in this respect see Ref. [7j. 

From the given A"-factors we derived the quark-mass dependence of helium nuclei us- 

n 

ing the techniques of Ref. [28]. Additionally we have presented a new response matrix of 
calculated primordial abundances to variations in nuclear binding energies and the scat- 
tering length of the important two-nucleon 1 S'o channel. Combining these we have derived 
a stringent limit on the quark-mass variation at the time of big bang nucleosynthesis of 
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Sm q /m q = 0.02 ± 0.04. 

In previous phenomenological studies the bounds derived from the deuteron and 4 He 
abundances are, e.g., 0.009 ± 0.19 and —0.005 ± 0.038, respectively, from Ref. [12j and 
—0.002 ±0.037 and 0.012 ±0.011, respectively, from Ref. [13]]. The uncertainties are signifi- 
cantly smaller compared to ours, since in these works no attempt was made to quanitify the 
theoretical uncertainty — in Ref. 12J it is stated that the uncertainty is expected to be of 
the order of a factor of 2. In this sense, although the uncertainty of our work seems larger, 
still the bound derived is more robust, since a careful uncertainty estimate was done. 



In Ref. |28j a range —1% < 5m q / m q < 0.7% is quoted for the quark mass variation allowed 
by BBN. In their calculation, the sensitivity is dominated by variation of the neutron lifetime, 
r, which strongly affects the 4 He abundances which are well-constrained observationally. 
We have shown that this calculation is based on a model-dependent assumption. Under the 
reasonable assumption that all elementary particle masses are proportional to v, the Higgs 
vacuum expectation value, the sensitivity of r to v is shown to be quite large (although 
only half that found by [28J]) and we obtain a limit \5v/v \ < 0.9%, which within this model 
translates into \5m q /m q \ < 0.9%. 



In Refs. 



12 



13] also bounds are derived from the 7 Li abundance that are different from 
zero, namely 5m q /m q = 0.016 ± 0.005 and Sm q /m q = 0.013 ± 0.002, respectively. Also 
here the uncertainties do not include the theoretical uncertainty of the input quantities. 
In our work we have not included the 7 Li abundance for several reasons. Firstly, the 7 Li 
abundance is very sensitive to variations in 7 Be and 7 Li binding energies, as well as A = 5 
resonances 13[ which have not yet been calculated using the methodology of this paper. In 
the future our study can be extended into this regime, as soon as systematic studies of the 
quark mass dependence of heavier nuclei are available, e.g., employing methods of nuclear 
lattice calculations — see Ref. 88|, |89J and references therein. In fact, the first results within 
that framework for 4 He, 8 Be and 12 C can be found in Secondly, the discrepancy between 
theory and observation is a factor of ~ 3 which can make nonlinear effects important. Lastly, 
the observational status and interpretation of the 'Spite plateau' of lithium abundances in 



Pop. Ill stars is still uncertain (see, e.g., Refs. 



83 



90 



91| and references therein). 



33 



Acknowledgments 



C. H. acknowledges the hospitality and the funding by UNSW where a part of this work 
was carried out. J. N. acknowledges funding by Fundacion Ramon Areces. This work is 
supported in part by the Spanish Research contracts FPA2007-29115-E/, FPA20 11-27853- 
C02-02, and FPA2008-05287-E/FPA "Quark Masses and Hadron Physics (From Quarks to 
life)", by the Australian Research Council, by the DFG and the NSFC through funds pro- 
vided to the Sino-German CRC 110 "Symmetries and the Emergence of Structure in QCD", 
by the HGF through funds provided to the Nuclear Astrophysics Virtual Institute NAVI 
(VH-VI-417), by the BMBF (06BN9006), by the EU via ERC project 259218 NuclearEFT 
and the EU FP7 HadronPhysics3 project. 



[1] J.-P. Uzan, Rev. Mod. Phys. 75, 403 (2003); Living Rev. Relativity 14, 2 (2011). 

[2] K. A. Olive, M. Pospelov, Y.-Z. Qian, A. Coc, M. Casse, and E. Vangioni-Flam, Phys. Rev. 

D 66, 045022 (2002) | hep-ph/0205269] . 
[3] R. L. Jaffe, A. Jenkins, and I. Kimchi, Phys. Rev. D 79, 065014 (2009) |arXiv:08Q9.1647l 

[hep-ph]]. 

[4] C. J. Hogan, Rev. Mod. Phys. 72, 1149 (2000) [astro-ph/990~9295l . 
[5] M. Livio, D. Hollowell, A. Weiss, and J. W. Truran, Nature 340, 281 (1989). 
[6] H. Oberhummer, A. Csoto, and H. Schlattl, Science 289, 88 (2000) [astro-ph/0007178| . 
[7] E. Epelbaum, H. Krebs, T. A. Lahde, D. Lee, and U.-G. Meifiner, larXiv:1212.4l8T1 [nucl-th] 
(2012). 

[8] J. K. Webb et ai, Phys. Rev. Lett. 107, 191101 (2011); J. A. King et al, Mon. Not. R. Astron. 
Soc. 422, 3370 (2012). 

[9] J. C. Berengut, V. V. Flambaum, J. A. King, S. J. Curran, and J. K. Webb, Phys. Rev. D 
83, 123506 (2011). 

[10] T. Dent, S. Stern, and C. Wetterich, Phys. Rev. D 76, 063513 (2007). 

[11] W. J. Marciano, Phys. Rev. Lett. 52, 489 (1984); X. Calmet and H. Fritzsch, Eur. Phys. J. C 
24, 639 (2002); P. Langacker, G. Segre, and M. J. Strassler, Phys. Lett. B 528, 121 (2002); 
C. Wetterich, J. Cosmol. Astropart. Phys. 10, 002 (2003); T. Dent and M. Fairbairn, Nucl. 

34 



Phys. B 653, 256 (2003). 
[12] V. V. Flambaum and R. B. Wiringa, Phys. Rev. C 76, 054002 (2007) |arXiv:0709.0077l [nucl- 
th]]. 

[13] J. C. Berengut, V. V. Flambaum, and V. F. Dmitriev, Phys. Lett. B 683, 114 (2010). 
[14] O. Civitarese, M. A. Moline, and M. E. Mosquera, Nucl. Phys. A 846, 157 (2010). 
[15] A. Coc, P. Descouvemont, K. A. Olive, J.-P. Uzan, and E. Vangioni, Phys. Rev. D 86, 043529 
(2012). 



[16 

[17; 

[18 

[19 
[20 

[21 

[22 
[23 
[24 

[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 



V. V. Flambaum and E. V. Shuryak, Phys. Rev. D 65, 103503 (2002) |hep-ph/0201303|; Phys. 



Rev. D 67, 083507 (2003) [hep^p h/0212403|. 

V. V. Flambaum, A. Holl, P. Jaikumar, C. D. Roberts, and S. V. Wright, Few Body Syst. 38, 
31 (2006) [nucl-th/0510"075] . 

A. Holl, P. Maris, C. D. Roberts, and S. V. Wright, Nucl. Phys. Proc. Suppl. 161, 87 (2006) 



|nucl-th/0512048| 



V. V. Flambaum, [physic s/0309107 (2003). 

V. V. Flambaum, D. B. Leinweber, A. W. Thomas, and R. D. Young, Phys. Rev. D 69, 115006 
(2004). 

M. E. Carrillo-Serrano, I. C. Cloet, K. Tsushima, A. W. Thomas, and I. R. Afnan, 
larXiv: 1208.30091 [nucl-th] (2012). 

H. Leutwyler, Annals Phys. 235, 165 (1994) | hep-ph/9311274] . 

P. C. Bruns and U.-G. Meifiner, Eur. Phys. J. C 40, 97 (2005) |hep-ph/0411223] . 

C. Hanhart, J. R. Pelaez, and G. Rios, Phys. Rev. Lett. 100, 152001 (2008) |arXiv:08Q1.287T1 

[hep-ph]]. 

J. R. Pelaez and G. Rios, Phys. Rev. D 82, 114002 (2010). 
J. Nebreda and J. R. Pelaez, Phys. Rev. D 81, 054035 (2010). 
J. Gasser and U.-G. Meifiner, Nucl. Phys. B 357, 90 (1991). 

P. F. Bedaque, T. Luu, L. Platter, Phys. Rev. C 83, 045803 (2011) [arXiv: 1012.38401 [nucl-th]]. 



E. Epelbaum, U.-G. Meifiner, W. Glockle, Nucl. Phys. A 714, 535 (2003) |nucl-th/0207089|. 



S. R. Beane and M. J. Savage, Nucl. Phys. A 713, 148 (2003). |hep-ph/0 206113|. 

E. Epelbaum, U. G.-Meifiner, W. Glockle, C. Elster, Phys. Rev. C 65, 044001 (2002). 

V. Bernard, N. Kaiser, and U.-G. Meifiner, Int. J. Mod. Phys. E 4, 193 (1995) 



|hep-ph/9501384|; Nucl. Phys. A 615, 483 (1997) |hep-ph/9611253| . 

35 



[33] J. Soto and J. Tarrus, Phys. Rev. C 85, 044001 (2012) [arXiv: 11 12.44261 [nucl-th]]. 
[34] J. Mondejar and J. Soto, Eur. Phys. J. A 32, 77 (2007) [nucl-th/0612"05l] . 
[35] J. F. Donoghue, Phys. Rev. C 74, 024002 (2006) l nucl-th/0603"016] . 
[36] J. F. Donoghue, Phys. Lett. B 643, 165 (2006) [nucl-th/0602"074"l . 
[37] S. Diirr et al, Science 322, 1224 (2008) |arXiv:0906.3599 [hep-lat]]. 

[38] A. Walker-Loud, "Baryons in/and Lattice QCD", plenary talk, Workshop on Chiral Dynamics 
2012, Jefferson Lab, Newport News, VA, USA (to appear in the proceedings). 

[39] M. Frink, U.-G. Meifiner, and I. Scheller, Eur. Phys. J. A 24, 395 (2005) | hep-lat/0501024] . 

[40] M. Albaladejo and J. A. Oiler, Phys. Rev. D 86, 034003 (2012) [arXiv: 1205.66061 [hep-ph]]. 

[41] J. M. Alarcon et al, arXiv: 1209.28701 [hep-ph] (2012). 

[42] J. Gasser, H. Leutwyler, and M. E. Sainio, Phys. Lett. B 253, 252 (1991). 

[43] C. Ditsche, M. Hoferichter, B. Kubis, and U.-G. Meifiner, JHEP 1206, 043 (2012) 
[arXiv: 1203.47581 [hep-ph]]; 

[44] M. Hoferichter, C. Ditsche, B. Kubis, and U.-G. Meifiner, JHEP 1206, 063 (2012) 
|arXiv:1204.625T1 [hep-ph]]. 

[45] G. S. Bali et al, Nucl. Phys. B 866, 1 (2013) |arXiv:1206.7034l [hep-lat]]. 

[46] A. S. Kronfeld, Ann. Rev. Nucl. Part. Sci. 62, 265 (2012) [arXiv: 1203. 12041 [hep-lat]]. 

[47] J. Bijnens, Prog. Part. Nucl. Phys. 58, 521 (2007). 

[48] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985); Nucl. Phys. B 250, 517 (1985); 

Nucl. Phys. B 250, 539 (1985). 
[49] G. Amoros, J. Bijnens, and P. Talavera, Nucl. Phys. B 602, 87 (2001). 
[50] G. Colangelo, J. Gasser, and H. Leutwyler, Nucl. Phys. B 603, 125 (2001). 
[51] A. Bazavov et al [MILC Collaboration], PoS CD 09, 007 (2009). 
[52] G. Colangelo et al, Eur. Phys. J. C 71, 1695 (2011). 

[53] C. W. Bernard et al, Phys. Rev. D 64, 054506 (2001); Ph. Boucaud et al [ETM Collabo- 
ration], Phys. Lett. B 650, 304 (2007); C. Allton et al [RBC and UKQCD Collaborations], 
Phys. Rev. D 76, 014504 (2007). 

[54] X. Feng, K. Jansen, and D. B. Renner, Phys. Rev. D 83, 094505 (2011) [ arXiv:1011.5288l 
[hep-lat]]. 

[55] S. Aoki et al, Phys. Rev. D 60, 114508 (1999). 

[56] C. R. Allton et al, Phys. Lett. B 628, 125 (2005); M. Gockeler et al [QCDSF Collaboration], 



36 



PoS (LATTICE 2008), 136 (2008) |arXiv:0810.5337l [hep-lat]]. 
[57] T. N. Truong, Phys. Rev. Lett. 61, 2526 (1988); Phys. Rev. Lett. 67, 2260 (1991); A. 
Dobado et al, Phys. Lett. B 235, 134 (1990); A. Dobado and J. R. Pelaez, Phys. Rev. D 
47, 4883 (1993); Phys. Rev. D 56, 3057 (1997). 

J. Nebreda, J. R. Pelaez, and G. Rios, Phys. Rev. D 83, 094011 (2011). 
E. Epelbaum, W. Glockle, U.-G. Meifiner, Eur. Phys. J. A 19, 125 (2004) |nucl-th/0304037 
J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984). 



[58 
[59 
[60 
[61 



[62 

[63 
[64 
[65 
[66 
[67 
[68 
[69 

[70 
[71 
[72 

[73 

[74 
[75 

[76 

[77 



R. G. Edwards et al. [LHPC Collaboration], Phys. Rev. Lett. 96, 052001 (2006) 



|hep-lat/0510062|. 

N. Fettes, Pion nucleon physics in chiral perturbation theory, PhD thesis, University Bonn, 
2000. 



V. Bernard, U.-G. Meifiner, Phys. Lett. B 639, 278 (2006) (hepdat/0605010]. 
V. Bernard, Prog. Part. Nucl. Phys. 60, 82 (2008). 

E. Epelbaum, W. Glockle, U.-G. Meifiner, Eur. Phys. J. A 19, 401 (2004) |nucl-th/0 308010 1. 

E. Epelbaum, W. Glockle, U.-G. Meifiner, Nucl. Phys. A747, 362 (2005) [ nucl-th/0405048] . 

R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

E. Braaten, H.-W. Hammer, Phys. Rev. Lett. 91, 102002 (2003) |nucl-th/030 3038|. 

E. Epelbaum, H.-W. Hammer, U.-G. Meifiner, A. Nogga, Eur. Phys. J. C 48, 169 (2006) 



|hep-ph/0602225|. 

E. Epelbaum, U.-G. Meifiner, and W. Glockle, (n^chth/0208040 (2002). 



S. R. Beane and M. J. Savage, Nucl. Phys. A 717, 91 (2003) |nucl-th/ 0208021 1. 

S. R. Beane, P. F. Bedaque, K. Orginos, and M. J. Savage, Phys. Rev. Lett. 97, 012001 (2006) 

[hep-lat/06020T0] . 

S. R. Beane, P. F. Bedaque, M. J. Savage, and U. van Kolck, Nucl. Phys. A 700, 377 (2002) 



|nucl-th/0104030|. 

J.-W. Chen et al, Phys. Rev. C 86, 054001 (2012) [arXiv: 1012.04531 [nucl-th]]. 

T. Yamazaki, K. -i. Ishikawa, Y. Kuramashi, and A. Ukawa, Phys. Rev. D 86, 074514 (2012) 

[arXiv: 1207.42771 [hep-lat]] . 

L. Kawano, "Lets go: Early universe II: Primordial nucleosynthesis the computer way," 
preprint FERMILAB-Pub-92/04-A (1992). 
C. Angulo et al, Nucl. Phys. A 656, 3 (1999). 

37 



[78] R. H. Cyburt, Phys. Rev. D 70, 023505 (2004). 

[79] S. Ando, R. H. Cyburt, S. W. Hong, and C. H. Hyun, Phys. Rev. C 74, 025809 (2006). 

[80] R. H. Cyburt and B. Davids, Phys. Rev. C 78, 064614 (2008). 

[81] E. Segre, Nuclei and Particles, (Benjamin-Cummings, Reading, MA, 1977). 

[82] D. Larson et al., Astrophys. J. Suppl. Ser. 192, 16 (2011). 

[83] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 010001 (2012). 

[84] R. H. Cyburt, B. D. Fields, and K. A. Olive, J. Cosmol. Astropart. Phys. 11, 012 (2008). 

[85] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 

[86] A. Walker-Loud, C. E. Carlson, and G. A. Miller, Phys. Rev. Lett. 108, 232301 (2012) 

|arXiv:1203.0254l [nucl-th]]. 
[87] A. Coc, N. J. Nunes, K. A. Olive, J.-P. Uzan, and E. Vangioni, Phys. Rev. D 76, 023511 

(2007). 

[88] E. Epelbaum, H. Krebs, D. Lee, and U.-G. Mei&ier, Phys. Rev. Lett. 106, 192501 (2011) 

|arXiv:1101.2547l [nucl-th]]. 
[89] E. Epelbaum, H. Krebs, T. Lahde, D. Lee, and U.-G. Meifiner, Phys. Rev. Lett. 109, 252501 

(2012) |arXiv:1208.1328l [nucl-th]]. 
[90] F. Iocco and M. Pato, Phys. Rev. Lett. 109, 021102 (2012). 

[91] K. A. Olive, P. Petitjean, E. Vangioni, and J. Silk, Mon. Not. R. Astron. Soc. 426, 1427 
(2012). 



38 



